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We analyze the behavior of an itinerant Fermi system near a charge nematic (n — 2) Pomeranchuk 
instability in terms of the Landau Fermi liquid (FL) theory. The main object of our study is the 
fully renormalized vertex function F n , related to the Landau interaction function. We derive T n 
for a model case of the long-range interaction in the nematic channel. Already within the Random 
Phase Approximation (RPA), the vertex is singular near the instability. The full vertex, obtained 
by resumming the ladder series composed of the RPA vertices, differs from the RPA result by a 
multiplicative renormalization factor Zr, related to the single-particle residue Z and effective mass 
renormalization m*/m. We employ the Pitaevski-Landau identities, which express the derivatives 
of the self-energy in terms of T n , to obtain and solve a set of coupled non- linear equations for Zr, 
Z, and 77i* /m. We show that near the transition the system enters a critical FL regime, where 
Zy ~ Z oc (1 + <? c ,2) 1,/2 and m*/m ~ 1/Z, where g C} 2 is the n = 2 charge Landau component which 
approaches —1 at the instability. We construct the Landau function of the critical FL and show 
that all but g c ,2 Landau components diverge at the critical point. We also show that in the critical 
regime the one-loop result for the self-energy 'S(K) oc J dPG(P)D(K — P) is asymptotically exact 
if one identifies the effective interaction D with the RPA form of T n . 

PACS numbers: 71.10.Hf, 71.27.+a 



I. INTRODUCTION 

The Fermi liquid (FL) theory states that low-energy 
excitations in a system of interacting fermions are rep- 
resented by fermionic quasiparticlcs which differ quanti- 
tatively but not qualitatively from free fermions.— 
Quantitative changes are encoded in the Landau func- 
tion, which is a tensor in the spin space and a function 
of the momenta of colliding quasiparticlcs in the orbital 
space. Angular harmonics of g in the charge and spin 
channels, g cn and g s ,n, determine renormalized values of 
various observables. 

The FL theory also allows for instabilities which occur 
either due to Cooper pairing or to symmetry-breaking 
deformations of the Fermi surface (FS). The latter are 
known as Pomeranchuk instabilities.— In the FL nota- 
tions, a Pomeranchuk instability occurs when one of the 
Landau harmonics approaches —1. 

Examples of Pomeranchuk instabilities include phase 
separation (g c .o = ~ 1) ; a ferromagnetic transition 
(g Si o = —1), at which the FSs of spin- up and spin-down 
fermions split apart, and nematic- type transitions in the 
chargeZ£&±^±^±i^^I and spin channels 
which lower the rotational symmetry of the FS. In mod- 
ern literature, the point in the parameter space where a 
Pomeranchuk instability occurs is called a quantum crit- 
ical point (QCP). 

As the system approaches a Pomeranchuk instability, 
a certain "FS susceptibility" , which measures "softness" 
of the Fermi surface with respect to deformations of par- 
ticular symmetry, diverges^^ In certain cases, e.g., for a 
ferromagnetic transition, a diverging FS susceptibility co- 
incides with the thermodynamic susceptibility whose di- 
vergence signals a phase transition. In general, however, 



the FS and thermodynamic susceptibilities are different. 
For example, any instability resulting in a deformation of 
the FS in an isotropic FL is not accompanied by a diver- 
gence of a thermodynamic susceptibility, which remains 
isotropic in the disordered phase. In those cases, a way 
to detect Pomeranchuk instabilities is to construct the 
correlation function of the incipient order parameter of 
certain symmetry* 7 - A susceptibility corresponding to this 
correlation function, Xa*,n*, diverges as 1/(1 + g a *,n*) at 
a QCP in the critical channel {a*,n*}, where a = c, s. 
Although it does not follow from any general principle, it 
is usually assumed that Landau components in all other 
channels arc not affected by a Pomeranchuk instability 
in the {a*,?i*} channel. 

Taken at face value, the last assumption im- 
plies that, near a Pomeranchuk QCP, the effective 
mass to* = m(l + g c ,i) remains finite. How- 
ever, Hertz-Millis-type effective theories of Pomeranchuk 

Qgp = ^9, 10, 12, 16, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34 m wmcn 

the original four-fcrmion interaction is replaced by the 
interaction of fermions with fluctuations of the incipient 
order parameter, predict a different behavior. In these 
theories, the effective mass diverges at a QCP in di- 
mensions D < 3: as | ln(l + g a *, n *)\ in D = 3 and as 
(1 + g a * ,n*) -1 ^ 2 m D = 2. As the FL theory is supposed 
to be valid at T = down to the very QCP, the diver- 
gence of to* implies that of g c ,i, in clear disagreement 
with the assumption that all but g a *,n* Landau compo- 
nents are not affected by the Pomeranchuk instability in 
a particular channel. 

As susceptibilities in the FL theory contain the fully 
renormalized mass, it then becomes an issue whether this 
singular mass renormalization should be included into the 
ordinary FL formula for the susceptibility in the critical 
{a*,n*} channel Xa* ,n» = (1 + 5c,i)/(l + g a *, n *) and also 
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whether non-critical channels can be "dragged" to criti- 
cality simply due to a divergence in the effective mass. 

In this communication, we re-visit this issue. We 
consider a Pomeranchuk QCP in the charge channel 
(a QCP in the spin channel deserves a separate con- 
sideration because of subtle issues related to spin con- 
servation, see Ref. [20t ). For definiteness, we focus on 
the charge nematic instability in the n = 2 chan- 
nel (Refs.0ii[T3[ll|[ll[ll[ll[ll[ll[I3); our conclusions, 
however, are also valid for all n > 1. We show that, 
for D < 3 the assumption of only one critical Landau 
component holds only at some "distance" from a QCP 
in the parameter space. In the immediate vicinity of a 
QCP, a FL of new type emerges. In what follows, we 
will refer to this new FL as to a "critical FL". The in- 
teractions at high energies, which were the cause of the 
Pomeranchuk instability in the first place, play the role 
of "bare" interactions for the critical FL. Accordingly, 
the enhanced nematic susceptibility Xc \ ^ 1/(1 + .9c, 2) 
is a bare susceptibility of the critical FL. The low-energy 
interactions, mediated by soft collective fluctuations in 
the n = 2 charge channel, lead to further renormaliza- 
tions of the FL parameters. These renormalizatons are 
encapsulated in the "critical "Landau function, g, which 
we show to have both charge and spin components with 
any n, even if the original, "high-energy" FL has only 
the n = 2 charge Landau component. 

Our key result is that all components of g diverge in 
the same way, i.e., as 1/(1 + g Ci -i) x ^ 2 , upon approaching a 
nematic QCP in D < 3. In particular, a divergence of g c ^\ 
implies that m*/m = l+g c ,i diverges at the QCP as well. 
At the same time, since divergences are the same for all 
g a ,m they cancel out in the expressions for the suscepti- 
bilities Xa.n <x (l+5c,i)/(l+5o,n)i which retain the same 
values Xa n as m the original FL. In particular, suscepti- 
bilities in channels different from the critical one remain 
finite at a QCP, while Xc,2 preserves its l/(l + g Ct 2) form. 
This means that the divergence of the effective mass near 
a Pomeranchuk QCP does not affect the behavior of any 
of the susceptibilities. 

To obtain these results, we derive diagrammatically an 
expression for the fully renormalized anti-symmctrized 
interaction vertex r n (k, p) between the particles with 
momenta k and p on the Fermi surface. This vertex 
is obtained from a more general vertex function in the 
limit of zero momentum transfer and vanishing frequency 
transfer and is the "input" parameter for the FL theory: 
the Landau function is proportional to T n 

g = 2 vZ 2 — f , (1.1) 

m 

where Z is the quasiparticle residue. 

We identify the most relevant part of T near a QCP 
and show that it describes an interaction mediated by 
soft collective bosonic fluctuations in the n = 2 charge 
channel. For small |k — p| w kp8, we find V (6) oc 
1/ [l + <7 Ci 2 + (akp9) 2 + ..] where dots stand for less rel- 
evant terms. The length scale a is the effective radius 



of the interaction U(q) in the d— wave charge channel, 
and the product afcp is a dimensionless parameter of our 
theory. The calculations are under control if ak-p ^> 1, 
which we assume to hold. For akp 3> 1, the system is 
the critical FL regime when (1 + g c ,2) < ^/(akp) 2 *C 1. 

Such a singular form of T was proposed earlier on phe- 
nomenological grounds 3 ^ and obtained within the Ran- 
dom Phase Approximation (RPA) for the Hubbard model 
near an antiferromagnetic instability— To obtain the full 
T n in our case, we first generalize the RPA result to 
the n = 2 charge Pomeranchuk instability. We show 
that the RPA-typc formula for the effective interaction 
f RPA (k, p) oc 1/(1 - U(k - p)n 2 (k, p)), where ILj(k, p) 
is the static polarization bubble in the n = 2 channel, is 
reproduced by summing up the diagrams which do not 
renormalize the bare interaction U. 

Next, we analyze the diagrams for T beyond RPA and 
show that in the critical FL regime the full T n differs from 
pQ.RPA by a cons t an t factor 1 ~ Z^ 1 . This relation 
results from re-summation of a particular non-RPA series 
of diagrams which includes renormalization of U into a 
full dynamic interaction. We show that all other non- 
RPA diagrams are relatively small in powers of 1/akp 
and can therefore be neglected. The existence of an extra 
factor Zr between T n and r fi,RPA is very important for 
our analysis - just the RPA form of T would produce 
nonsensical results in the FL description. 

Another input parameter for a FL theory is the quasi- 
particle residue Z. To obtain Z, we use the exact 
Pitaevski-Landau relation s 1 ! 37 which express Z in terms 
of T n (k, u)k', p, u) p ), where now p and k are not necessarily 
at the FS and Wk.p arc finite. We extend the previous cal- 
culation of r to fermions away from the FS and obtain 
Z as a function of (1 + <? c ,2) 1 ^ 2 with akp as a parameter. 
This function is rather complex but reduces to a simple 
form Z ~ (akp)(l + gc^) 1 / 2 in the critical FL regime. Us- 
ing this form of Z and full T n , we construct the Landau 
function of the critical FL and show that all its compo- 
nents diverge in a way discussed above. We also show 
that the effective mass m* diverges at but not before 
a QCP. In this respect, our results do not support the 
conjecture 3 ^ that the effective mass may diverge before 
the system reaches a QCP. 

We also discuss the relation between the exact formula 
for the self-energy S(/c, u>) to linear order in lu and k — k-p, 
obtained from the Pitaevski-Landau relations, and a one- 
loop formula for £ due to an exchange by soft bosonic 
collective excitations. We show that the one-loop for- 
mula is asymptotically exact in the critical FL regime if 
the effective interaction is replaced by ZT Q , i.e., by the 
RPA form of the effective vertex p^ RPA . The correc- 
tions to the one- loop formula are small in 1/ak-p and in 
M/o>fl, where wfl oc (1 + g c ,2) 3 ' 2 is the upper bound- 
ary of the FL behavior. We emphasize that the one-loop 
approximation is valid only for the linear-in-w term in 
the self-energy. The next, to 2 In a; term has contributions 
from all orders even if the theory is extended to a large 
number of fcrmionic flavors 7Vi^ Outside the FL regime, 
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i.e, for \lo\ > wfl, the self-energy scales as w 2 / 3 . Some 
vertex corrections in this regime are small in 1/N (Refs. 
25 |Hl32|33:), while others remain O(l) even for large N 

(Ref. iiir 

Regarding the full form of the self-energy near a QCP, 
we show that, for 1 + g c ,2 <C 1, the self-energy is "local", 
in the sense that it depends primarily on to but not on 
(-k = v-p(k — kp). The prcfactor for the e k term in the 
self-energy scales as 1/afcp and is, therefore, small. 

A Pomeranchuk instability in the d— wave charge chan- 
nel was introduced in the context of the RG analysis of 
potential instabilities of a 2D Hubbard model^ Shortly 
thereafter, Oganesyan et al. analyzed in detail a d— wave 
charge instability in isotropic systems.— The subject at- 
tracted substantial interest^^'l 2 -^'^^^ both from 
the theoretical perspective and also due to potential rel- 
evance to cuprates^ and ruthenites.— There are some 
subtle differences between a d— wave Pomeranchuk in- 
stability in lattice systems ^ 10 ' 11 and an n = 2 Pomer- 
anchuk instability in isotropic systemS ) 9 ' 12 ] 13 i 14 i 15 i 16 but 
our conclusion is that the physics does not change qual- 
itatively between isotropic and lattice cases (see below). 

Properties of the Fermi liquid near a charge Pomer- 
anchuk instability were studied by Rosch and WolfleJ^ 
They did not consider the Landau function near a ne- 
matic QCP, but obtained the effective mass in terms of 
the critical parameter l + g c ,2 using the one- loop approxi- 
mation for the self-energy with an RPA form of the effec- 
tive vertex. Our results agree with Ref. [l5| in that the ef- 
fective mass can be obtained within the one-loop approx- 
imation with the RPA vertex, but our dependence of the 
mass on l + g c ,2 is different from that in Ref. [H. The dis- 
agreement originates from the difference in the forms of 
the static polarization bubble for dressed fermions, which 
we discuss in Sec. IIIIB II 

The structure of the paper is as follows. In Sec. HH we 



review a general procedure of constructing the FL ver- 
tices. Section IIHI is devoted to the critical FL theory. 
In Sec. IIII Al we introduce a model for the nematic 
Pomeranchuk instability in the charge channel. In Sec. 
IIII Bl we obtain the FL vertex T n near this instability 
within the RPA approximation. In Sec. IIII CI we go be- 
yond the RPA level and obtain the full vertex T n . In 
Sec. IIVI we obtain the quasiparticle residue Z, the ver- 
tex renormalization Zp. and effective mass m* from the 
exact Pitaevski-Landau relations. In the same Section, 
we also analyze the crossover between ordinary and crit- 
ical FLs. In Sec|Vj we obtain the Landau function in 
the critical FL using the full vertex found in Sec. IIII CI 
and show that all components of this Landau function di- 
verge at a QCP. This is the main result of our paper. In 
Sec lVIl we discuss the relation between the one-loop and 
exact Pitavcskii-Landau forms of the self-energy near the 
FS. We present our concluding remarks in Sec. IVIII A 
number of technical issues are discussed in Appendices 
EHB andO 

II. DIAGRAMMATIC DESCRIPTION OF AN 
ORDINARY FERMI LIQUID 



We begin with a brief overview of the diagrammatic 
description of an ordinary FL. The FL theory describes 
effects of the interactions between fermions confined to 
a near vicinity of the FS. The interactions which in- 
volve fermions away from the Fermi surface are ab- 
sorbed into the Landau function g a p,-yS- This function 
is related to an exact antisymmetrized vertex T^j 1 $ via 
5a/3,7(5 = 2^Z 2 (m*/m)r2 / 3 j7 ,5 {v = m/2n is the density of 
states in 2D). The vertex T^g s is defined for particles 
at the FS in the limit of zero momentum transfer q and 
vanishing energy transfer CI (Ref. [H) 



r2 ft7a = Jim Q [r(k, 0; p, 0|k - q, -fi; p + q, fl)S ai 6 0S - T(k, 0; p, 0|p + q, fi; k - q, -fl)5 aS S 0J ] \ k=P =k F , (2.1) 

I 



where T(k, w^; p, w p |k, w^; p, w p ) is an exact nonsym- 
metrized vertex. The momentum transfer is equal to q in 
the first term of Eq. (|2.ip and to k p q in the second 
one. 

For a generic ratio of q and f2, diagrams for the ver- 
tex r(k, Uk', p, w p |k, uJk', P, Wp) [and its exchange coun- 
terpart r(k, <x>fc; p, u)p\p, L) p ] k, u>k)] can be separated into 
two groups depending on whether they contain propaga- 
tors of "soft" particle-hole pairs 

P(q,fi,fc) = -vj de k J ^G(k + q,^ + ft)G(k,^) 

(2.2) 



with q — ► and r2 — ^ 0, where e k = v-p(k — k-p). The 
meaning of V is the propagator of a particle-hole excita- 
tion with momentum q and energy f2 formed by fermions 
moving in the direction of k = k/fc. 

For vanishingly small q and f2, the dynamic part of the 
bubble V is determined by fermions on the FS; therefore, 
diagrams with soft bubbles are attributes of the FL the- 
ory. For these diagrams, the order of limits q — > and 
O — * matters because V = in the limit of q/fl — > 
and is finite (equal to v) in the limit Q/q — > 0. The di- 
agrams without soft bubbles generally involve fermions 
with high energies, of order E-p. For these diagrams, the 
order of limits q — > and f2 — > is irrelevant. 
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Because r^. 7 a is the vertex in the limit q/Q, — > 0, 
it does not include soft particle-hole bubbles (V = in 
this limit) and, in that sense, it is a high-energy property 
playing a role of the bare vertex in the FL theory. The 
second-order diagram for r9g s are shown in Fig. [Da). 
The wavy line in these diagrams is the bare, static in- 
teraction potential V(k, p;k,p) which may depend not 
only on the momentum transfer k — k but also on the 
incoming momenta themselves. 

On the other hand, physical observables, e.g., the spe- 
cific heat or the nematic charge susceptibility, contain 



contributions from soft bubbles because the observables 
are affected by elastic collisions between particles right 
on the FS, i.e., by processes with ft — and finite q. The 
corresponding antisymmetrized vertex is called 5 . 
In Fig. [if)), we show the second-order diagrams for 

r a/3 ;7 5- Thc ful1 vcrtcx r */3; 7 5 is obtained from r^. 7(5 by 
summing up an infinite series of ladder diagrams shown 
in Fig. [T]e), which leads to a familiar integral equation 1 



J 



/d6 



(2.3) 



where a shorthand Kp,Pp denotes the "four- 
momentum" on the FS, i.e., Kp = {kp,Wfc = 0}, 
kp = fcpfc, and 9 is the angle between k and k'. 

The full vertex r^. 5 is related to the scattering am- 
plitude fa/3-fS in thc same way as T^g 5 is related to 
g a f3;~f8, i-e., via f a/ 3 n s = 2vZ 2 (m* /m)T g ag . nS . The re- 
lation between / and g takes a particularly simple form 
when expressed via partial components in the charge and 
spin channels, /; and gf. 



fa,l 



9a,l 



9a,l 



(2.4) 



Using these relations, one obtains FL formulas for ob- 
servables in terms of g a ,i- 



III. CRITICAL FERMI-LIQUID THEORY 
A. Nematic charge fluctuations 




K ^ P 



P K 



We follow early wor kWWV 7 and consider a ne- 
matic charge instability described by the model Hamilto- 
nian with the interaction in the d— wave charge channel 



H = Yj ^(q) rf k VU q /2,a C p-q/2,/3 C P+q/2,/3 C k-q/2,, 
k,p,q 



Here 



<ik = \/2cos(20k) 



(3.1) 
(3.2) 



FIG. 1: Fermi-liquid vertices to second order in the bare in- 



teraction (wavy line), a) Diagrams for V 



for U 



n 

a/3;7<5 • 



b) Diagrams 



between V 



n 

a/3;7(5 



Diagrammatic representation of the relation 
and T 9 a g 5 , Eq. (|2.3p . A bubble composed 
of dotted lines represents the particle-hole propagator V [see 
Eq. (|2.2|| ] in the limit Q/q — > 0, where V reduces to a con- 
stant. 



space, i.e., 



is the "d-wave" formfactor and 0k is the angle between 
k and an arbitrarily chosen x-axis. Hamiltonian (|3.ip 
describes the interaction between nematic fluctuations of 
the electron density p d (q) = Ek rf kcj c+q /2, Q c k-q/2,Q- 

To keep the treatment under control, we assume that 
the interaction U(q) is sufficiently long-ranged in real 



U(q) = U P(qa), 



(3.3) 



where the function P(x), subject to P(0) = 1, is a de- 
creasing function of its argument, and thc effective inter- 
action radius a is much larger than the inverse average 
distance between particles, i.e. kpa 1. We also assume 
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that U(q) is analytic for small q: 
U(q) = U (l - (qa) 2 



(3.4) 



Equation (|3.1[) is a reduced version of a more general 
interaction between quadrupolar fluctuations of the elec- 
tron density cf a QjjC^ ai where Qij = <5yV 2 — Idjdj, and 
i,j = x, y (Rcfs. I9lll7lh Such an interaction can be de- 
coupled into a longitudinal part, which is the same as 
in Eq. (|3.ip . and a transverse part with the momentum- 
dependent factor d^dp, where dk = v / 2sin(20k)- The two 
terms contribute separately to the fermionic self-energy 
and r° and can be treated independently of each other. 
To simplify the presentation, we consider only the lon- 
gitudinal part of the quadrupolc interaction, given by 
Eq. (|3.ip . The effect of the transverse part on thermo- 
dynamic properties of a FL has recently been considered 
by Zacharias, Garst, and Wolfleii. 

The nematic charge susceptibility of free fermions is 
defined as 



xS(q^) = 2n d (q,fi) 

= -2^d k G(k + q/2,w fe 



(3.5) 

n/2)G(k-q/2,Wfc-n/2). 



We will be interested in the long wavelength and low 
frequency d-density fluctuations with q <C kp and 
|f2|/upg -C f. In this regime, xfl f° r f ree 2D fermions 
with quadratic spectrum £k = k 2 /2m is 

,2 



xS(q,fi)= 21/(1 



^-2cos 2 (20 q )H 
2k£ H v F q 



(.3.6) 



Notice that Landau damping of nematic fluctuations is 
anisotropic whereas the dispersion of static xfl with q 
(absent for n = density fluctuations for q < 2fcp) is 
isotropic. 

To first order in the interaction, * s gi ven by the 

first two diagrams in Fig. [TJ a): 

r 2/3; 7 5( k :P) = UodfrdpSaydps - f/(k - p)d1 +p 5 a s6 g ~ / . 

(3.7) 

According to our definition of dk, when both k and p are 
on the FS, 



cos(2</>k) cos(2(^ p ) 



1 — sin(20k) sin(2</> p ) 

(3.8) 



where dots stand for non-cos(2</>k) terms which we ne- 
glect. Using Eq. jSl]) and the SU{2) identity 



(3.9) 



we separate Eq. ()3.7[> into the charge and spin channels 



r? = d k d. 



(3.10) 



[/b--tf(k-p) 



r'j = --d k dpU(k-p). 



a) 

X c ^ — ^^^^ + + ^^^^ — 



K K 






r 4 



FIG. 2: RPA series for the nematic susceptibility [a)] and 
Solid circles represent d-wave vertices dk and 
wavy lines represent the forward-scattering part of the in- 
teraction, Uo- The bubbles are with zero vanishing external 
momentum and zero frequency -the same bubbles as in Fig. [T] 
6,c. Since, by assumption ak^ <C I, only the diagrams that 
contain Uo are included. 



At this stage, the interaction is static and mostly in 
the d-wave channel. Hence the effective mass equals to 
the bare one and the quasiparticle residue Z = 1. The 
Landau function is then obtained by simply multiplying 
Eq. (|3.10|) by 2v. Since, by assumption, J7(k — p) is 
peaked at k = p, the charge d— wave Landau harmonic, 
g c ,2-, is much larger than the spin d-wave harmonic g s ,2- 
Indeed, for ak-p ^> 1, 



9c,2 = 2z/ 
9s, 2 = 2i/ 



-^d k d p T c ps 2isU | 1 



2?r 2tt 
d</> k d^>p 



c 



2tt 2tt 



■d k d p r s 



-2i/l7 — — < 5 C 2 , 
afcp 

(3.11) 



where c = (3/47r) Jg 00 dxP(x). To leading order in 1/kpa, 
we then have g Cj 2 = 2vU$, g Sj 2 = 0. We see that g c ^. 
approaches —1 when Uq approaches — l/(2i>). 

With only g c ^ being non-zero, the nematic spin sus- 
ceptibility Xs,2 = X°l(l + 5c,i)/(l + 9s,2) = xf,l retains 
its bare value, while the charge susceptibility is 



(o) 1 +9c,i 

Xc,2 = xA- 



(0) 
Ac, 2 



1 + 9c, 2 1 + 9c,2 



(3.12) 



diverges when g c ^ tends to —1. This corresponds to the 
"conventional" scenario of a Pomeranchuk QCP, where 
only one of the Landau components approaches —1 but 
the effective mass does not diverge. We will see be- 
low that such a behavior does not survive in a non- 
perturbative theory with fully renormalized vertices. 

Before moving to the analysis of vertex renormaliza- 
tion, we note that Eq. (|3. 12[1 and the FL relation (|2.4[) 
can be reproduced diagrammatically by calculating di- 
rectly the charge susceptibility by calculating r*g In 
both cases we need to include only those diagrams that 
contain n<i(q, fi) at vanishing momentum and zero fre- 
quency and, thus, do not contribute to rcnormalization 
°f Summing up the ladder series for Xc,2 in Fig. [5] 
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FIG. 3: RPA diagrams for F 



The difference with the 



series for T q a p s in Fig. [2^ is in that the polarization bubbles 
are evaluated at finite external momentum k p. 



i.e., to retain only diagrams with a maximum number 
of polarization bubbles Il^q, 0) at any given order, now 
with finite q = k p. While all such diagrams contain 
the interaction only in the form U(k — p) , non-RPA dia- 
grams involve integrals of U over intermediate momenta. 
Since each of these integrals contributes a small factor of 
1/akp [see Eq. (|3.11[) ]. non-RPA diagrams can be safely 
neglected. Summing up the RPA series and neglecting 
the first order Uq term, which is irrelevant near a QCP, 
we obtain 



we obtain 



Xc, 2 (q^ 0,0 = 0) = 



(0) 

Ac, 2 



y (0) 
Ac, 2 



1 + 2vU 1 + g c 



which coincides with Eq. (|3.12p . The ladder series for 
r a/3; 7 5 in Fig. [2] 6 yields 



U 



1 + 2isU 



(3.14) 



Using the relation between r*^. 7i5 and the scattering am- 
plitude /a/?,7,5 = 2vT q a/3rfS , we find 



fc 



2U u 
l + 2U is' 



(3.15) 



Comparing Eqs. (|3.11|) and (|3. 15[> . we see that the FL 
relation 



/< 



c,2 



5c,2 



l+.9c 



(3.16) 



T n ^ A (k,p 



C/(k- 



a/3; 7 5 \*,VJ~ 2 1 + 2c/(k _ p)nd(k _ p Q) S aS S 0y . 

(3.17) 

Notice a "wrong" combination of spin indices (S a s5/3 7 



(3.13) 

instead of 8 ai 8ps)- 



Thermodynamic parameters of a FL (including the ef- 
fective mass) are determined by r^^J^(k, P) f° r ^ ne P ar ~ 
tides on the FS and at zero frequencies. However, for 
the consideration in the next sections, we will also need 
to know r^ R ^ away from the FS, i.e., for momenta 
different from k-p and for non-zero frequencies of incom- 
ing and outgoing fcrmions. A simplification is that, for 
1 + .9c.2 "C 1, relevant momenta are still close to hp, while 
the relevant values of |k — p| and uj^, u> v are small. In this 
situation, F^^^ is still given by an expression similar 



to Eq. (j3~17l) 



T a kjS (k,w fc ;p,u;p) 



1 



d k d p 8 a s8f3 7 U(k - p) 



2 1 + 2U(k - p)II d (k - p, w fc - u p ) 

(3.18) 

but where now H c i(k — p, tjft — Up) is a dynamic bubble. 



is indeed reproduced. 



The polarization bubble 



B. Random Phase Approximation for nematic 
instability 

At the next step, we consider renormalization of 
r^. 7< 5(k, p). An assumption of the long-range interac- 
tion in the nematic channel allows one to use the RPA, 



The d-wave bubble of free fermions for q <C kp and 
"C vpq: is given by Eq. (|3.6|) . Substituting this form 
into Eq. (|3.18|) . replacing U(k — p) by Uq in the numer- 
ator, and neglecting the q 2 /2kp term compared to the 
(aq) 2 one which comes from the expansion of U(k — p) 
in Eq. (|3.4|) . we obtain 



Pn P ^l+ 5ci2 + |k-pp a 2 + 2cos2(2^ q )i^^ 
= d^ 5 ai 8, s + a ai -a ps 

8^ 1 + g c 2 + | k - p|2 a 2 + 2 C0s2(2^ q )l|^ 

I 

where, as before, </> q is the angle between q and an arbi- This form of r^ R ^(k, p) is only valid, however, at 

trary chosen x-axis. some distance away from a QCP, where fermions behave 

as nearly free quasiparticles. Near a QCP, the quasipar- 



1 1 



7 



ticle mass and residue differ from the free-fermion values 
and II ( j(q, 0) is to be treated as a fully renormalized bub- 
ble, which we label as II*, (q, f2). 

Renormalization of the polarization bubble is a sub- 
tle issue and we pause here to discuss it in some detail. 
In general, the polarization bubble contains two types of 
renormalizations: the self-energy insertions, which trans- 
form the bare Green's functions into the renormalized 
ones, and vertex corrections. We will show in Sec. IIII C 31 
that the vertex corrections are small in the quasistatic 
limit (fi <C vpq); thus, in this limit, the bubble can be 
computed as a convolution of two renormalized Green's 
functions. Still, one has to be careful even with this com- 



putation because the quasiparticlc residue Z and the ef- 
fective mass m* depend on the energy they are measured 
at: near the FS, they approach the renormalized values 
of Z < 1 and m* > m, while at higher energies they 
approach the free values Z = 1 and m* = m. It turns 
out that renormalizations of the three terms in the ex- 
pansion of the free bubble [Eq. (|3.6[) ] are determined by 
different energies. First, we consider the constant term 
IT2(q — > 0,Q = 0). For free fermions, it coincides with 
the density of states. A product of two Green's functions 
with close arguments can be represented as a sum of two 
parts: coherent and incoherent (Ref. |l|): 



G(k + q, w + Q)G(k, u) = 2nZ 2 J F ' q - 6(w)8(e k ) + *(k, «), (3.20) 

vj, • q — ill 

I 



where = vp(m/m*)k. A contribution of the coherent 
part to n^(q — > 0,17 = 0) is equal to the renormalized 
density of states vZ 2 (rn* /m). However, it would be in- 
correct to replace IT^q — > 0, Q = 0) only by the coherent 
contribution because there is also a contribution from the 
incoherent part, <&(k, w), which cannot be evaluated ex- 
plicitly. A way avoid this complication is is to integrate 
over the momentum first. This integral comes from high 
energies, of order of the ultraviolet cutoff of the theory, 
A. At these energies, one can approximate G by its free- 
fermion form G _1 (k, oj) = to — e k . Integrating over in 
between —A and A first and then over the frequency, we 
find that n^(q, 0) is given by the bare density of states: 

U* d (q - 0, A = 0) = -v I p [ A de k 1 

J 2n J-A (iuj - e k ) 

[duo A , 
= " / 2— T2 = 3 -21 

Both energies (tu and e^) in this integral are of order A, 
which justifies the replacement of the Green's function by 
its free-fermion form. This result is known in the theory 
of electron-phonon interaction^ renormalization of the 
phonon spectrum by particle-hole excitations is the same 
as if the fermions were free. 

Similarly, the q 2 term in Eq. (|3.6[) is also a high-energy 
contribution. If the electron spectrum remains quadratic 



up to A, this term is the same as for free fermions; other- 
wise, hp in the prcfactor of this term is replaced by some 
momentum of order of cither kp or of the reciprocal lat- 
tice spacing. 

In distinction to the static part of II^, the dynamic 
part of the bubble- the Landau damping term-comes 
from low-energy fermions by virtue of energy conserva- 
tion. Calculating this part using only the coherent term 
in Eq. (|3.20[) . we see that the Landau damping term is 
multiplied by an overall factor of (Zm* /m) 2 . 

The full II*,(q, £1) is then given by 

n3(q,^)=,f 1 -^-2cos 2 (2^ q )f^l) 2 B), 

(3.22) 

where b ~ kp. 

The q 2 - and Landau damping terms in Eq. (|3.22|) agree 
with their counterparts in the expression for the polariza- 
tion bubble in Ref. [l^. The first, constant term differs, 
however, from that in Ref. ITU where only the contribu- 
tion from coherent fermions was included and, as a result, 
the first term had an extra Z factor. 

Substituting Eq. (|3~2"2"]l into Eq. (f3~T%)) and neglecting 
again the q 2 /b 2 term in the bubble compared to the {aq) 2 
from the interaction, we obtain the dynamic effective in- 
teraction as 



r O.RPA/i \ 

W k, Wfc ;p,. p) - 4, 1+5c)2 + |k _ p|2a2 + 2cos2(2 ^ q)m ^ 

= d k dp5 a s60 7 T(k - p,u k - w p ). (3.23) 
I 



1 d^dpSasSpj 



For k and p near the FS, |k — p| = 2kp sm(0/2) where |k — p| « kpO. Notice that by virtue of Eq. (|3.9[) . the ef- 
9 = Z(k, p). Because k and p are nearly parallel, 
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fective interaction has both charge and spin components 
even though we neglected the g s ^ component of the bare 
interaction. 

A cautionary note: Eq. (|3.23|) is obtained within the 
RPA, which does not mix channels with different angu- 
lar momenta. The RPA is valid, strictly speaking, only 
for a long-range interaction which is not a very realistic 
assumption. For a generic interaction, there is no proof 
that the effective interaction does not change quantita- 
tively beyond the RPA, e.g., it is possible that non-RPA 
renormalizations destroy a simple pole structure of Eq. 
(|3.23p . In Appendix [Al we discuss the diagrammatic se- 
ries for Y n beyond RPA and show that the effective in- 
teraction does change beyond the RPA. 



is true, however, only if the interaction is static. Due to 
Landau damping, the dressed interaction has a singular- 
ity in the complex plane, i,e., a branch cut in the Matsub- 
ara formalism. In this situation, the frequency integral of 
the product G 2 (q, fi) times the dressed interaction has an 
extra contribution from the Landau damping singularity. 
To account for this effect, we have to re-consider dia- 
grams with soft particle-hole bubbles and replace each 
static interaction, U(k — p), by a dressed one, i.e., by 



Y 



k, Uk',p,u} p ) given by Eq. (|3.23[) . Since the RPA 
interaction already contains all insertions of particle-hole 
bubbles, the remaining diagrams form the ladder series 
shown in Fig. gj In Sees. IIIIC2I and IIIIC3I we evalu- 
ate this ladder series and show that corrections to ladder 
diagrams are small. 



C. Diagrammatics of the critical FL theory 

1. Interaction via dynamic collective fluctuations 

We now show that the RPA expression for & 
[Eq. (|3.23p ] is not the full result, even if the control pa- 
rameter akp is large4^ The reason is that the diagram- 
matic series for ^ in the bare interaction, considered 
in the previous Section, does not include diagrams with 
particle-hole bubbles at exactly zero momentum and van- 
ishingly small frequency. The argument is that such dia- 
grams vanish due to particle-number conservation. This 



2. Renormalization ofT^g s beyond 



RPA 



A building block for ladder diagrams is the product 
of two Green's functions with the same momentum and 
frequency and the dynamic interaction (|3.23p . We will 
analyze explicitly the first few terms in the ladder series 
and then sum up infinite series of diagrams. Consider 
first the second-order diagram in Fig.|4]a). Summing over 
the internal spin indices and neglecting the momentum 
transfer q compared to kp in the d-wave formfactors, we 
obtain 



J 



rg(k,p)=4^ 



g/^/^(k + q ,0)f( qi O)f(k-p + q ,f2) ) 



(3.24) 



r 



where |k| = |p| = kp, = ^(q, k), and f is defined 
in the second line of Eq. (|3.23|) . Let us integrate over 
first. Both the Green's functions and f depend on 
0. The (^-dependence of Y is in two places: first, in the 
anisotropic Landau damping term, which depends on 
as cos 2 [2 (0t + </>)], where 0k = Z(k, x), and, second, in 
the magnitude of the momentum |k — p + q|. We will see, 
however, that the integral over is dominated by narrow 
regions near ±7r/2, where q is almost perpendicular to 
k. Therefore, the prcfactor of the Landau damping term 
can be replaced by cos 2 (20k). Expanding next near 
±7r/2 as = ±7r/2 =p with |0 <C 1, and using the fact 
that near a QCP the relevant values of |k — p| « kpO 
(determined by the inverse correlation length) are small, 
we find that |k — p + q| 2 w (kp9 ± q) does not depend 
on to leading order. It is thus only G 2 (k + q, 51) that 
depends on 0. The integral over can then be evaluated; 
one has to be careful though to resolve an ambiguity at 



Q = 0, which results in a delta-function term: 



J-t 2tt y v F q 



mi 



^2l3/2 



m [n 2 + {v F qy 

(3.25) 

The delta-function term will give the leading contribution 
to the vertex. It is reproduced if the integral over is 
performed as 



2tt J-oo 2ir fan 



v F qq 



5(ft) 
v ¥ q ' 



(3.26) 



The second term in Eq. (|3.25p is reproduced if one keeps 
track of the actual limits of the integral over 0. 

Substituting Eq. (|3.25[) into Eq. (|3.24|1 and integrating 
over Q and q, we find, as advertised, that the leading 
contribution to r^.^-(k, p) comes from the S— functional 
term, while the contribution from the second term in 
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a) 

K 



K 



P K 



K 



P K 



P K 



K 



a/3;7<5 



beyond 



(|3.25p is proportional to 1 + g c .2, which is small near a 
QCP. The delta- function of Q eliminates the frequency 
dependence of the RPA vertices, and, integrating over q, 
we obtain for the second-order ladder diagram 



^n,{2} 

a/3;7<5 



(k, P ) 



dud, 



4 



— 5a6$, 



Z 2 m*X 
2m 



Av 



4(1 + g c ,2) + (01*0)' 



,(3.27) 



where 



A = 



2ak F ^J\ + 9c,2 

We will see below that Z 2 m*\/m « 1 in the critical 
FL regime. This implies ^ is of the same order as 
pn.RPA anc ^ ^hus, ladder renormalizations cannot be ne- 
glected. We also note that the parameter A will play the 
key role in our further analysis. In particular, it will be 
shown that the critical FL regime is defined by a con- 
dition A <C 1 rather than by 1 + g c ,2 <C 1. Because 
A contains a small factor of l/ak F <C 1, the former is 
stronger than the latter. 

We now return to the ladder series composed from the 
G 2 (q, f2) blocks. Repeating the same procedure for the 
nth term of the series, we obtain 



r n ' {n} fk n) 



d 2(n-l) 
dkd p — 

' Z 2 m*\ 



d^ 



2m 



2 

n-l 



where 



4;/ 



n 2 (l+g c , 2 ) + (ak F 9y 



fM(6»), (3.29) 



(3.30) 



Comparing Eq. (|3.30p to Eq. (|3.23[) . we see that the struc- 
ture of the RPA interaction (the n = 1) term is not re- 
produced at higher orders: first, the dependence on the 



angle 9 between k and p is different and, second, there 
are additional d factors arising at each order of the per- 
turbation theory. However, the angular dependence can 
be scaled back to that of the n = 1 term by replac- 
ing 9 — ► n9. This means that the angular harmonics 
[ <MH™}(0)cos(£0)/27r of f{"> arc the same as those of 
r, as long as t and n are not too large: In <C A. Since 
A <C 1 in the critical FL regime, this condition is always 
satisfied for not too large I and n. In other words, f^™} 
can be re-written as the sum of two terms: the first one 
is the same as for n = 1 and the second one vanishes on 
angular integration 



FIG. 4: a) Ladder series for the full vertex F 
RPA. Each hatched block represents T^p.^g at the RPA level, 
Eq. (|3.23|l . &)) One-loop vertex correction to the ladder dia- where 
gram. 



r<">(0) =T(0) + rW(0), (3.31) 



R {n} {9) = (n-l). 



(ak F 9) 2 - n(l + g ca ) 



[n 2 (l 4- g c , 2 ) + (ak F 9) 2 } [1 4- g c ,2 + (ak F 9) 2 } 

(3.32) 

and / d9R^ = 0. Since the obscrvablcs are determined 
by harmonics of r n , we will neglect the remainder term 
below, i.e, replace f^{9) by f (6). 
Still, the subsequent terms in the series contain ad- 
ditional factors d^™ -1 "* + dp . These terms account 
for renormalization of a smooth, regular variation of P 
along the FS. Such a variation, inherent to a model 
with anisotropic, d-wavc, cos(2</>) interaction, implies 
that both the quasiparticle residue and effective mass 
are not uniform along the FS. The c?-channel vertex can 
be obtained by extracting contributions proportional to 
(3.28) dh.d p at any order, neglecting all other terms, and sum- 
ming up the series. On the other hand, one can average 
the vertex over the FS, i.e., restore Galilean invariancc 
broken by neglecting the sin(2</>) interaction. The two ap- 
proaches lead to physically equivalent results which dif- 
fer only by numerical prefactors. For brevity, we only 
discuss here the approach based on averaging, which is 
particularly appealing because it allows one to use the 
technique developed for isotropic FL systems to obtain 
to* /m and Z. Results of alternative approaches are pre- 
sented in Appendix [Bl We emphasize that, although no 
important physics is lost in cither of the approaches, all 
of them arc still approximate. We will discuss one subtle 
issue associated with this approximation in Appendix. [U1 
We now apply the averaging procedure. Since relevant 
k and p are almost parallel to each other, we will set 
k = p in the d factors and use an identity 



(d 



2ns (2n-l)!! 



(3.33) 



Substituting this expression into L^'^ 71 ^ and summing up 
the series, we find that the ladder series for T^p s result 

in overall renormalization of r^/^f/^ 



1 -pO.RPA 

1 ct/3;7<5 



1 1 



Z T Avl + g C:2 + [ak F 9) 2 ' 

(3.34) 
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where 



— = V 
z T ^ 



n=l 



Z 2 m* 
2m 



(2n-l)!! 



z 1 , 



X [1 + Jl 



Z-i 



-X 



.(3.35) 



Since Z F < 1, ladder renormalization enhances the inter- 
action compared to the RPA result. 



and introducing polar coordinates x = r cos iJj, y 
rsinip. We then obtain 

fm*\ 1/2 1 
r A q^ 0,0 = ) = -Z — —x 

\m J ak F 



where 



— dr dtp-„ 
4tt ./o Jo > cos W + a + tan V 



sin tp cos ip 



" = (l + ffc,2)- 



, (3.39) 



(3.40) 



Integrating over r, we obtain 



3. Vertex corrections 

We now show that vertex corrections to ladder series 
are small and thus Eq. (|3.34[) is a complete result for our 
model. The lowest order vertex correction to T n is shown 
in Fig. 0J>. When calculating the ladder diagram in the 
previous Section, we saw that the interaction vertices are 
effectively static (O = 0), while typical bosonic momenta 
Q ~ \/l + Qc.il a are small near a QCP. To estimate the 
vertex correction, we can then simply calculate a three- 
leg vertex, shown in Fig. \Sp, at zero external frequency 
and finite but small external momentum q. With these 
simplifications, the three-leg vertex reduces to 



r A (q^ 0,0 = 0) 



dfl' d 2 q' 
~2T(2^r) 



xG(k 



rG(k + q',0') 

q,0')f(q',0'). 

(3.36) 



Z fm*\ 1/2 1 



r A (q^ 0,0 = 0) = -- — -j-H<*), (3-41) 
8 \ m J ak F 



where 



For x < 1, 



v/2 , , sin^ 
dip =. 
o vx + tan?/' 



(3.42) 



0.847 -0.927x- 



(3.43) 



In the next Section, we show that a ~ 1 away from a 
QCP and a <C 1 near a QCP. In both cases, J-(x) ~ 1. 
Therefore, the vertex correction can be estimated as 



(3.44) 



We neglected the spin and <i-wave factors as well as the 
anisotropy of the Landau damping term, all of which give 
only overall numerical prefactors. Since T(q',n') is now 
isotropic, the angular integration in Eq. (|3.36[) can be 
performed first. Because the momenta in two Green's 
functions differ by small yet finite q, the angular integral 
of G(k + q', 0')G(k + q' + q, O') contains only a regular 
but no 5(0') term, i.e., the result is given by the second 
term in Eq. (|3.25p . Hence 



which is at most of order 1/afcp <C 1. 

Equation (|3.44p can be re-written via the character- 
istic momentum q or characteristic frequency O from 
Eq. (j335)l as 



r A ( q -> o, o = o) ~ qo /k F ~ n /E F 



F! 



(3.45) 



r A ( q ^ o,o-o) 



d9! dq'q' 



2?r 2tt 
xf(q',0'). 



where E F = k F v F . When written in this form, it is clear 
that r A (q — ► 0,0 = 0) is an effective Migdal (adia- 
batic) parameter, i.e., the ratio of characteristic energies 
of bosons (O ) and fermions (Ep), of our theory^ For 
ak F 3> 1, this ratio is small both away and near a QCP. 
[(O') 2 + (v F q') 2 ] 3 ^ 2 The Migdal parameter occurs here in same way as for 
the case of electrons interacting with optical phonons. 
Namely, if the interaction in Eq. (|3.36|) is replaced by a 



10' 



(3.37) 



We notice immediately that the remaining double inte- 
gral is finite even right at a QCP, where l+g c .2 = 0. This 
is to be contrasted with ladder diagrams of the previous 
section which diverge near a QCP. An explicit calculation 
can be carried out by rescaling the variables as q' = xqo 
and O' = Oq?/, where 



Zy/m*/m/a, O = q a v F , 



(3.38) 



propagator of optical phonons Oq/(0 2 + O 2 ) (where now 
Oo E F is the optical phonon frequency), the result for 
the vertex is the same (up to a number). The reason 
for the smallness of the vertex in both cases is a small 
phase space defined by q and O . Although the vertex 
correction is also small for acoustic phonons, the reason 
is different in this case: a typical momentum of acous- 
tic phonon is of the order of the inverse lattice spacing 
(~ k F ), while the electron energy at the same momentum 
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is of the order of Ep, i.e., much larger than the phonon 
energy - 

Two cautionary notes, both related to the fact that 
the integral in Eq. (|3.41[) is determined by 0' ~ v F q' ~ 
(E F /(ak F )Z^/m/m*. First, we used a quasistatic, 
Q'/v F q' form of the Landau damping term, which is 
not, strictly speaking, justified in this regime. How- 
ever, because q' is not much smaller than Q'/v F ei- 
ther, the result for Ta is still correct up to a numer- 
ical prefactor. Second, and more important, we as- 
sumed a renormalized FL form of the Green's function 
G(k, <jj) = Z/[u> — v F (k — k F )]. Meanwhile, we will see 
later that (i) a FL behavior extends only up to fre- 
quencies of order ujfl ~ { v f/o){^ + 9c, 2) 3 ^ 2 , and (ii) 
ft ~ Z^J^(v F /a) ~ {vp/a){akpfl 2 {\ + g c , 2 ) 3/4 - 
Thus, typical Q' ~ flo <C lofl are outside the FL regime 
for small enough 1 + g c> 2,. However, the non-FL effects 
are not important here either because the three-leg vertex 
outside the FL regime is determined by very high energies 
(~ vp/a), where fermions behave as almost free quasi- 
particles with Z,m*/m w 1 (Ref. HH). The result is that 
Eq. (|3.45p still holds, but the relevant fermionic energy 
is Ep while the relevant bosonic energy is 11q ^ Ep/cikp 7 
and the vertex correction is simply given by 



r A (q^ 0,0 = 0) - l/akp. 



(3.46) 



Certainly, the same result for the vertex is obtained if one 
substitutes the free Green's function and bare vertex into 



Eq. (|3.44p from the very beginning. We also emphasize 
that the vertex correction was found to be small in the 
limit Cl/q — ► 0. In the opposite limit, q/Vt — > 0, vertex 
corrections are not small and must be included into the 
theory to ensure that the Ward identities are satisfied. 



Complete form of the FL vertex 



To summarize this Section, we conclude that, in a crit- 
ical FL, the full vertex s differs from the RPA re- 
sult by a constant renormalization factor 1/Zp given by 
Eq. (|3.35j) . Hence a complete result for the static d— wave 
vertex, averaged over the FS, and for k « p reads 



a/3;7<5 



(k, P ) 



1 



01 



AvZp 1 + g c a + {akpOf 



(3.47) 



To obtain an explicit form of r™g s , we need to know Z 
and m* jra. This is what we discuss in the next Section. 
We show there that, near the QCP, the leading terms in 
Z and m*/m are determined by the static vertex, i.e., 
there is no need to invoke dynamics of bosonic modes. 
The sublcading terms, however, require the knowledge of 
the dynamic vertex. Following the same steps that led us 
to Eq. (|3.47[) . one can show that, for akpd 3> */l + g c ^, 
the dynamic vertex is given by 



J 



W * ;P -" P> = ^l + 9 , J + (aW + C(^) 2 ^' 



(3.48) 



where C ~ 1. In a general case, the dynamic vertex is a 
scaling function of the following form 



J 



1 Srv/iS, 



AvZp 1 + g C: 2 ' 



ak\ 



'6 ( Zm*\ a|wfe — ujp\ 
~g c ^ , \ m J v F {l + g c , 2 ) 3/2 



(3.49) 



where Q(x, 0) = 1/(1 + x 2 ). As an explicit form of Q 
will not be required for the estimates of the sublcading 
terms in Z and m*/m, we will use below the approximate 
Eq. (PH81) with C = 1. 



IV. PITAEVSKI-LANDAU RELATIONS AND 
SELF-CONSISTENT EQUATIONS FOR Z AND 

m* /m 

A. Pitaevski-Landau identities for the derivatives 
of the self-energy 



The renormalized mass can be found in two ways. If 
the Landau function is known, m* /m can be determined 
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from its g c ^\ component. Alternatively, one can compute 
the self-energy in a given microscopic model and extract 
mass renormalization from the renormalized Green's 
function. In general, the quasiparticle residue Z can be 
found only from a microscopic calculation because renor- 
malization of Z is determined by states away from the FS, 
which are not described by the FL theory. Renormaliza- 
tions of Z and m* /m near QCPs have been considered 

j n 9 . 10. 12. 16. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34 y j & l 00 p_ w i se 

expansions in the effective interaction. Near a QCP, how- 
ever, typical energies involved in renormalizations van- 
ishe in the inverse proportion to the divergent correlation 
length. This allows one to determine both m*/m and Z 
within the FL framework. To do so, we will use the 
Pitaevski-Landau identities which relate the derivatives 
of the seli-energyiS! to T n from (JS3Z1). 

We will need three of the Pitaevski-Landau identities 

i-E/^(^p)[ G ^)] 9 |f 



■A 

d 3 p 

n (2tt) 3 ' 
(4.1) 



Here, as before, a shorthand K F denotes the "four- 
momentum" on the FS, i.e. K F = {kpk,ujk = 0}, while 
the objects [G 2 (P)] f2 and [G 2 (P)] <] represent the prod- 
uct G(P)G(P + Q) taken at vanishing Q in the limits of 
q/il — > and 0,/q — > 0, respectively. Notice that the in- 
ternal momentum P is generally not at the FS. Also, we 
switched here from the Matsubara formalism, employed 
in the previous sections, to the causal one. To simplify 
notations, we use the same symbols for Matsubara and 
real frequencies in all cases when it does not lead to a 
confusion. 



The first two relations are essentially the Ward identi- 
ties following from particle-number conservation (gauge 
invariance), while the third one is special for Galilean- 
invariant system. We remind that Galilean invariance 
was restored in our problem by averaging over the FS. 
The vertices T n and T q are related by Eq. IpOj) , 



Equations (|4.ip and (|2.3p determine a Taylor expan- 
sion of the self-energy S(fc,cj) to first order in u> and 
4 = v F (k-k F ): 



J 



E(#) = iE 





-(w-e*) / T%, af) (K F ,P) [G 2 (P)} 



d 3 P 



k I "p£7 



fi (2tt) 3 Z 
I 



T%^{K F ,P) [G 2 (P)]-[G 2 (P)] i 



p • k F d 3 P ' 

kj (2tt) 3 . 

m 



where the FS term T,(K F ) is absorbed into a shift of the valid for any order of integration over the fermionic mo- 
chcmical potential. Using an identity mentum and frequency, we can rewrite Eq. (|4.2p as 



[G 2 (K)]-[G 2 (K)] n = - 2 -^P^S(c 



kp 



(4.3) 



E ( JQ = (w - e k ) 







Jr2^(K F ,p) [c 2 (p)] n ^ 



Zm 
4^2 



E 





a0,a0 



(0) cos ede. 



(4.4) 



Equivalcntly, 



where 



E(A-) = (w - e k ) 



1 

1 

Z 



(_k_ 

Z 



m 

1 

m* 



(4.5) 



T% iaft (K F ,P) [G 2 (P)] ^ 



d 3 P 
(2tt) 3 



(4.6) 
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1 

m* 



1 

m 



47T 2 



(4.7) 



This is an implicit equation for m* /m. 



and 8 is the angle between k and p when both momenta 
are on the FS. 

The integral in Eq. (|4.7[) involves a full static vertex of 
the interaction between the particles on the FS. In our 
case, this vertex is given by Eq. (|3.34|) with lo^ = uj p . 
Using this equation, we find 



7 2 

— Y 

47r 2 t-r 1 



r% ta0 (6) C o S OdO = A- A 



and 



m 



= 1 



A Z 2 m* 
4 Zy m 



(4.8) 



Next, wc consider the Z factor. In contrast to the 
effective mass, the Z factor is, in general, determined by 
the dynamic vertex. However, we will see shortly that 
the leading term in Z still comes from the static vertex. 
To see this, we substitute the dynamic interaction from 
Eq. (|3.48|) into Eq. (|4.6[) and re-label the variables as 
(4.9) q = p — kp and fi = uj p , upon which Eq. (|4.6[) reduces to 



1 iZ 2 

— = 1 

Z AZ T v 



d 2 qdfl 



(2tt) 3 (ft + Q- 



r 



ft-o i + ffc , 2 + (a(7)2 _(^l) 2 ! M' 



(4.10) 



This integral is similar to that in Eq. (|3.24[) for the 
second term in the ladder series for s , except for 
now we have one rather than two interaction vertices. 
The result of the angular integration is the same as in 
Eq. (|3~25|) . The leading contribution to l/Z comes from 
the ^-function term in Eq. (|3.25[) . which renders the ver- 
tex static. The contribution from the second, dynamic 
term Eq. (|3.25p is proportional to the three-leg vertex, 
estimated in Sec. IIII C 31 except for an extra factor of 
1/Zy because the Z factor contains a full rather than 
RPA vertex. Since the three-leg vertex is small, so is the 
dynamic correction to \/Z. Collecting all terms, we find 



1 A Z TTi* m* . . 

7=1 + T ^ r A , (4.11) 

Z 4 Zy m m 



where Ta ~ Ta ~ 1/akp <C 1 is the dynamic correction 
due to a three-leg vertex. 



B. Solutions for m*/m, Z, and Zr near a QCP 



Equations l|4.9p . (|4.11|) . and (|3.35|) form a closed set 
from which one obtains m*/m, Z and Zy as functions of 
the coupling constant A, defined by Eq. (|3.28p . Neglect- 
ing A in Eq. (|4.11[) and comparing it to Eq. (|4.9[) and, 
we see that 



1. 



m 



(4.12) 



Using this relation, we can re- write Eqs. (|4. 1 1[) and (|3.35|) 

as 



1 4 / 1 -Z 
'Zy'^^Z 1 \z 



\Z 



1 _ 1 

~Z~ + 2 VI - AZ (1 + VI - A^) 



(4.13) 
(4.14) 



Introducing a new variable x = \Jl — \Z, we re- write Eq. 
(|4.14p as a cubic equation for x: 

a; 3 + a; 2 + (2A - l).x = 1. (4.15) 

The only real solution of this equation is 



l £ _2(A-2/3)_ 

3 L 1 ' 



(4.16) 



1/3 



where L = [9X + 8 + 3V3AV8A 2 - 13A + 16) . Once 

x is known, the three parameters of the FL theory 
m*/m, Z, and Zy are also known. For A <C 1, x = 
l-A/2 + 0(A 3 ); for A > 1, x= 1/2A+ 1/4A 2 + C(l/A 4 ). 
Using these asymptotics, we immediately find that at 
weak coupling 

771 

— L«i = l+V*; Z\ X<<1 = l-A/4; Z r \ x<<1 = 1-3A/4, 

(4.17) 

while at strong coupling 



7/7 



^L» 1 = A + 1 /4A; ;Z| a>>1 = 1/A-1/4A 3 ; 



m 
Zy 



1/4A+ 1/4A 2 



(4.18) 
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Notice that Z F ps Z/A in the strong coupling limit. 

We also emphasize that our calculations show that the 
efective mass to* diverges right at the critical point but 
not earlier. In this respect, our results agree with Ref. l29l 
but disagree with Ref. H3, where it was conjectured that 
the effective mass may diverge at a topological transition 
which pre-empts a QCP. Such a behavior would follow 
from our formalism only if, for some reason, renormaliza- 
tion of Z and Z F occurrs in such a way that Cz = Z 2 jZ F 
were an invariant. According to Eq. (|4.9p . the effective 
mass to*/to = (1 — AC^/4) would then diverge at finite 
A = i/Cz- However, the solution of the full set of equa- 
tions shows that the invariant is Z/Zt, in which case the 
divergence of to* is only possible at A = oo. 



Frequency and momentum dependences of the 
self-energy 



We now return to the self-energy given by Eq. 
In the previous section, we found that the relation 
Zm* I'm = 1 holds as long as the vertex correction A 
is small. Substituting this relation into Eq. (|4.5p . we find 
that £ is "local"- it depends on to but not on e^: 



The efe-dependence of S is determined by the vertex 
correction A. Substituting a full form of of 1/Z from 
Eq. (|4~TTjl into Eq. (|4~5]) . we find 



E(0,e fc ,) =e k 1- 



TO* Z 



T A e k ~ (4.20) 



We see that the k-dependent part of the self-energy re- 
mains small near a QCP. 



CRITICAL FL THEORY: 
PARAMETERS 



LANDAU 



£(#) ps (m*/m - l)w « Xlo. 



Explicit solutions for to*, Z, and Zp, obtained in the 
previous Section, complete our task of obtaining the FL 
vertex: one just needs to substitute these solutions into 
Eq. (|3.47[) . Having found the full vertex, we can now 
construct the FL theory of the critical region, i.e., re- 
late harmonics of 5 to observables. We emphasize 
again that T^p. s plays the role of an effective interaction 
between quasiparticles of an ordinary FL with "bare" ne- 
matic susceptibility Xc,2 oc 1/(1 +g c ,2) already enhanced 
by the interactions. We label "bare" susceptibilities as 
Xa,n and the ones renormalized by ^ a p.^s as Xa,n- 

Ap- 



(4.19) plying the standard FL phenomenology, we obtain 



Introducing the correlation length £ = a/jl + g c .ii 
we can re-write the mass renormalization coefficient as 
A = £/fcpa 2 . Notice that £ is the bare correlation length 
which enters the RPA formula for the vertex. Rosch and 
Woffle^, obtained the same expression for to* /to in terms 
of £ but, for reasons displayed in Sec. IIIIB 11 their £ con- 
tains an additional factor of \j\fZ compared to ours. 



= 1 + 9c\ 

TO 

1 +9cl _ _ 1 + gel , . 

Xc,n — Xc,n 1 _ j Xs,n — Xs,n _ , ^O.IJ 

1 + 5c,n 1 + g S ,n 

where <? Ci „ and g s ^ n are charge and spin components of 
the Landau function for the critical FL 



9ap n m = 2vZ 2 —T% n5 {6) 



Z 2 m* 



1 



4Z r m 1 + g c ,2 + {ak F 6) 2 1 1 p ' 1 

I 



g c ,2 + (ak F 9) 2 



(6 al 5f3S + Part ■ Vps) 

^21 



In 2D, the harmonics of g(0) are given by g n = numerator cancels with that of the effective "g-factors" 
§ g(8) cos(n9)d9/ (2ir). Integrating over 9, we find that 
the first n < A harmonics of g both in the charge and 
spin channels diverge concurrently with the effective mass 
upon approaching QCP as 



_ 1 + g c , 1 _ 5m _ 

Xa,n — Xa,n -. _ ~ Xa,,n _ ~ Xa,' 
r ~r ga,n ga,n 



(5.4) 



A. 



(5.3) 



all 



where A is given by p.28p . At a QCP, where A — > 
harmonics of g diverge. 

We also see, however, the susceptibilities retain their 
bare values despite the divergence of the Landau compo- 
nents. Indeed, because all Landau components diverge in 
the same way, renormalization of the effective mass in the 



In particular, the nematic susceptibility remains equal 
to Xc,2 oc 1/(1 + g c ,2), i.e., it is not affected by mass 
renormalization. 

One can also introduce an "actual" Landau function, 
g* , which describes a combined effect of rcnormalizations 
in the ordinary and critical FLs and includes contribu- 
tions from both the regular and collective mode parts of 
f^g. £. The nematic susceptibility can be written equiv- 
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FIG. 5: Top: schematic behavior of the charge (c) and spin 
(s) components of the critical Landau function g as a function 
of g c ,2- Bottom: same for the actual Landau function g* . 



The nematic component of the actual Landau function 
approaches —1 at a QCP, where g c .2 = — 1. We show the 
behavior of g a ,n and g* n in Fig. 

Note in passing that, although Xc,2 diverges as 1/(1 + 
g c ,2) m both the ordinary and critical FL regimes, this 
divergence has different origins in the two regimes. In 
an ordinary FL (where m*/m ~ 1), the divergence of 
\c,2 is entirely due to that of the d— wave charge "g- 
factor" 1/(1 + 17^2), while in the critical FL the divergence 
is equally "shared" between the effective mass and the 
"g-factor", each contributing a factor of 1/(1 + gc^) 1 ^ 2 - 
The analogous "sharing" holds for the spin susceptibility 
near a ferromagnetic instability. 



VI. EQUIVALENCE OF THE ONE-LOOP AND 
EXACT RESULTS FOR THE SELF-ENERGY 

In effective low-energy theories of QCPs, e.g., in the 
spin-fcrmion model, the effective interaction is D(Q) = 
9fbX{Q)i where x(Q) is the susceptibility of the diver- 
gent order parameter and (?fb is the fermion-boson cou- 
pling. The self-energy is obtained via a loop-wise ex- 
pansion in D, which is usually truncated at the one-loop 
order (cf. Fig. [TJa) 



alently either in terms of either g or g*: 

_ _ 1 + gc.i _ Xc°> } 1 + ffc.r _ {0} 1 + g*,l 
Xc2 ~ Xc,2 1 + ^ 2 - l + 1 + a - Xc.2 l+g * 2 

(5.5) 

All other susceptibilities Xa,n with {a, n} different from 
{c, 2} can be written as 



v _ v {o} 1 + g^i _ fo> 1 + g*,i 

Aa,n A a ,n i , - A a ,n i i _* 

1 T ga,n 1 T y a ,n 



(5.6) 



Comparing the expressions for x i n terms of g and g* , we 
read off the components of g* as 



9*a,n =9a,n + 9a,n{^ + ga 1 n). 



(5.7) 



For all partial components different from the n = 2 
charge one, the regular contribution <? a n is absent, and 
9a n = 9a.n- The divergence of g a ,n implies that all com- 
ponents of the actual Landau function g* different from 
the nematic one diverge at a Pomeranchuk instability. 
The n = 2 component behaves as 

9*c, 2 = (l + ffc,a)(l + ffc,a) - 1 



Sil(A') = J D{K-P)G{P) 



d 3 P 
(2tt) 3 



(6.1) 



Such a formula is used in the Eliashberg and FLEX the- 
ories (with the bare or full G(P), respectivel y 44 ! 45 ) . 

On the other hand, we showed in the previous Section 
that the linear in to and Ck parts of the E can be found 
from the Pitaevski-Landau identities using a fully renor- 
malized vertex. Now we can ask the following question: 
how does the one-loop result for E correspond to that 
obtained from the Pitaevski-Landau identities? In this 
Section, we show that Eq. (|6.ip is asymptotically exact 
in the critical FL regime (with corrections small in 1 /akp 
and w/wfl, if the effective interaction D(Q) is identified 
with ZT n (more accurately, with Y^n Z^^p a p)- 

To prove this assertion, we substitute J^p ZT^p a p for 
D into Eq. (jrTTj) . use Eq. (pTTT)) for T n , replace G by its 
coherent part, and take the limit of small uj and e^. The 
effective interaction can be decomposed into static and 
dynamic parts. The static part gives an "anomalous" 
efe term in Eil, which comes from the immediate vicin- 
ity of the FS. The dynamic part gives a regular ui — 
term, which comes from the entire phase space. Explic- 



itly, E 1L (/v) 



16 



saw 
Kim 



Z 2 m 



Air 2 ^ 



V TO* / 



H)E / r ^(A'F,P) [g 2 (p) 



d 3 p 

n (2tt) 3 ' 



(6.2a) 
(6.2b) 



r 



We sec immediately that the w terms in Eqs. (|4.4p and same if 
(I6.2al6.2b|) are the same, whereas the terms are the 



Z(Z- l)m 
4^2 



E 



a/3, a/3 



(0)cos 



E / 

/3 



r 



r/3 p 

(W) [G 2 (P)] ^3 



0. 



(6.3) 



Expressing the integrals of the vertices in Eq. (|6.3[) via 
to* and i? using Eqs. (|4.6[) and (|4.7p . we see that Eq. 
(|6.3p reduces to an identity 



Vto 



1 



-(Z-l) 



1 



m * 



(6.4) 



We thus see that both to and e*, terms in the one-loop self- 
energy coincide with the exact expressions obtained using 
the Pitaevski-Landau identities, if ZT n is identified with 
the effective interaction. This equivalence indeed holds 
only as long as to* /to and Z are energy independent and 
only a coherent part of G(Q) is relevant, which is the 
case for the critical FL regime of a QCP. 

In Appendix [C] we show explicitly how the self-energy 
is reproduced by the diagrammatic cxapnsion for r . 



VII. CONCLUSIONS 

In this paper, we analyzed properties of the Fermi 
liquid near a quantum phase transition, using a simple 
model of the nematic n=2 charge Pomcranchuk insta- 
bility as an example. Our main result is that, near a 
phase transition, the system enters into a new critical 
FL regime, in which all spin Landau components and 
all charge components with n ^ 2 increase and eventu- 
ally diverge at the critical point. This behavior is the 
consequence of a singular momentum dependence of the 
Landau function in a critical FL. Therefore, a common 
assumption that all Landau components, other than the 
one corresponding to the critical channel (g c ,2 for the ne- 
matic charge instability), are featureless near a transition 
is incorrect in D < 3. The divergence of the Landau com- 
ponents, including the one controlling renormalization of 
the effective mass, has no consequences for susceptibili- 
ties channels because the divergent effective mass cancels 
out with the divergent effective "g-factor" . 



To prove these statements, we derived the Landau 
function for a critical Fermi liquid, related to the ver- 
tex T n via g = 2i>Z 2 (m*/m)r n , where v is the density 
of states, Z is the quasiparticle residue, to* is the effec- 
tive mass, and F n is defined in the limit of zero momen- 
tum transfer and vanishingly small energy transfer. Our 
starting point is a model of 2D fcrmions with a d— wave 
interaction, which we assumed to be of sufficiently long 
range, such that ak-p ^S> 1, where a is the effective ra- 
dius of the interaction. We computed the Fermi-liquid 
vertex T a in two stages. First, we considered only those 
diagrams for T n that do not contain soft particle-hole 
bubbles with zero momentum and vanishing frequency. 
For ak-p 3> 1, diagrams of that type can be summed up 
in the RPA. The RPA vertex r n ' RPA contains a com- 
ponent which diverges at the critical point and can be 
interpreted as arising from the exchange by soft collec- 
tive excitations in the n = 2 charge channel. Next, we 
included the diagrams with soft particle-hole bubbles. In 
an ordinary FL, the full interaction is approximated by 
its static term, and such diagrams vanish due to particle- 
number conservation. In a critical FL, however, the effec- 
tive interaction is dynamic, and the diagrams with soft 
particle-hole bubbles are finite. We showed that the rele- 
vant diagrams of that kind form a ladder series. Summing 
up this series, we found that non-RPA diagrams with soft 
bubbles renormalizc the RPA vertex by a constant factor, 
i.e., the full vertex T n is equal to r n ' RPA /Zr, where Zp is 
a function of to* /to, and of the dimcnsionless coupling 
constant A = l/(2afcpwT + g c .2)- Vertex corrections to 
ladder diagrams were shown to be small as 1/akp. Hav- 
ing found r n , we employed the Pitaevski-Landau identi- 
ties for the derivatives of the self-energy in terms of T n 
to obtain coupled equations for Z and to* /to. This gave 
us a closed set of equations for Zp, Z, and m* /in with A 
as a parameter. Solutions of these equations allows one 
to follow the evolution from an ordinary FL behavior at 
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A < 1, where m*/m, Z and Zy are all close to 1, to a crit- 
ical FL behavior at A ^> 1, where m*/m = 1/Z « A and 
Zy ~ Z . We also showed that the self-energy of a critical 
FL is "local", i.e., S(fc, cj) « and that m* diverges 

only at A = oo rather than at finite A. The latter result 
implies that no preemptive transitions occur before the 
ncmatic one is reached. 

Finally, we showed that the exact self-energy (as given 
by the Pitaevski-Landau identities) is asymptotically 
close to the one-loop expression oc J G(Q)D{K — 

Q) if one identifies the effective interaction D with ZY = 
pn.RPA } iave ghown explicitly that in the limit of 

lu, €k —* 0, corrections to the one-loop result are small 
in powers of 1/akp and w/wfl, where ujfl is an upper 
bondary for the FL behavior. 
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it can be interpreted as the exchange by soft collective 
excitations. 

The issue of validity of the RPA is not specific to a 
charge nematic QCP. To avoid unnecessary complications 
associated with the d— wave factors in the vertices, we 
consider here a simpler case of an instability in the n = 

k k ( p Kryp 
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APPENDIX A: DIAGRAMMATIC SERIES FOR 
r" BEYOND RPA 



In the main part of the text, we adopted a model case 
of the long-range interaction in the the d-wave channel 
to justify the RPA. Such a model, although a necessary 
pre-requisite for an RPA-likc treatment in any channel, 
is hardly realistic and one is naturally led to wonder as 
to what extent the main features of the effective inter- 
action, Eq. (|3.23p , survive beyond the RPA level. In 
this Appendix, we address this issue by evaluating T n in 
a direct order-by-order perturbation theory for a short- 
range, Hubbard- like interaction (U =const). Notice that 
this issue is different from renormalization of the RPA 
interaction via ladder series in the dynamic interaction, 
considered in Scc. lIII Cl of the main text. Here, we discuss 
the full r without contributions from the soft partcle- 
hole bubbles, which is the input for the calculations in 
Sec. IIII CI The issue is whether this full T n has the same 
structure as the RPA result, Eq. (I3.23[) . i.e., whether 



FIG. 6: a) Third-order diagrams for r^g.„g that contain only 
the particle-particle and particle-hole bubbles at combina- 
tions of external momenta K ± P. b) Examples of the third- 
order diagrams that contain convolutions of particle-particle 
and particle-hole bubbles with Green's functions over the in- 
ternal momenta. 



channel (charge instability for U < or spin instability 
for U > 0). 

Diagrams for T n to second order are shown in Fig. [1] 
a). Explicitly, 



18 



T%. nS (K 7 P) = --(!- Un ph {K -P) + Un pp {K + P)) SasSfr 

+ - (1 + UU ph (K - P) + UU PP (K + P)) a aS ap y , (Al) 

I 



where IL ph (K) = - f(d 3 L/(2ir) 3 )G L G L +K and 
n pp (K) = - f(d 3 L/(2n) 3 )G L G- L+K are particle- 
hole and particle-particle bubbles. Notice that, to 
this order, all contributions to r^g.,^ contain either 
Tlph(K — P) or Ylp P (K + P) but no bubbles with other 
momenta. 

The new physics emerges at the third order of the inter- 



action. The third-order diagrams can be divided into two 
classes. Diagrams from the first class, shown in Fig. |5|)a 
still contain either H p h(K — P) or Ii pp {K + P). Dia- 
grams from the second class contain convolutions of Tl p h 
and Upp. 

Adding the contributions of the third order diagrams 
from the first class to Eq. (|A1[) . we obtain 



U 



T%. lS {K, P) = --(i-U (U ph (K -P)- n pp (K + P)) + U 2 (U 2 ph (K - P) + U 2 pp (K + P))) 6 aS S Pl 



u 



1 + U (U ph (K -P) + UTlppiK + P)) + U 2 (U 2 ph {K -P) + U 2 (K + P))) a aS ■ S 01 



(A2) 



This result for .~${K, P) can be re- written as an ex- pansion of the following formula 



(K,P) 



U 
~ ~2 

u 

+ 2 



1 



1 - UTlpp(K 
1 



1 



P) 1 + UUphi-P) 
1 



1 - UU pp (K 



p) i - un ph (K - p] 
I 



- 1 
- 1 



<?a8 ■ Vf3~f 



(A3) 



Although the last formula is valid only to order U 3 , it is 
very likely that higher-order diagrams for T n of the same 
type, i.e., containing only TL p h(K— P) or H pp (K + P) , are 
described by this expression. We see from Eq. (|A3[) that 
the structure of the vertex is virtually the same as in the 
RPA in a sense that F^p-iS contains three separate ge- 
ometric series, describing exchange by charge, spin, and 
pairing fluctuations. The three separate contributions 
to T^p-ys diverge near a corresponding QCP, e.g., near 
UU ph (0) = -1 for a charge QCP and UU ph {0) = 1 for a 
spin QCP (the former can occur if pairing instability is 
suppressed by e.g., magnetic field). We emphasize that, 
at this level, the particle-particle and particle-hole chan- 



nels entirely decouple, i.e., the interaction in the pairing 
channel does not affect the structure of the effective in- 
teraction mediated by soft collective excitations in the 
particle-hole channel. 

The situation is changed by diagrams from the second 
class. There are 24 topologically distinct third-order di- 
agrams which contains momentum integrals integrals of 
the polarization bubbles, and we refrain from presenting 
all of them. As an example, two of the 24 diagrams are 
shown in Fig. [H] b). We computed analytically all 24 di- 
agrams. There are numerous cancellations between the 
diagrams, and the final result is rather compact: 



a/3;7<5 



(K,P) = r^(K,P)-U 



d 3 L 
(2^)3 



[H p h(L) (2Gl+kGl+p — Gl+kGp-l) + U pp (L)Gl-pGl-k] S a sSpj 
f d 3 L 

J (27r) 3 ^ph^^L+KGp-L + H pp (L)Gl-pGl-k] &aS ■ 0>y, (A4) 
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where T^.^g is given by Eq. 



We see from Eq. (|A4[) that there are cross terms 
which involve both particle-hole and particle-particle 
bubbles, i.e., the terms with H p h(L)GL+KGp-L and 
H pp (L)Gl-pGl-k)- The presence of such terms im- 
plies that the particle-particle and particle-hole channels 
do couple beginning from the third order in U . This 
coupling should affect the interaction mediated by near- 
critical charge- and spin-density fluctuations. 

Note that T contains a term which involves only 
particle-hole bubbles [(ILj/j {X)Gx+kGx+p\- To under- 
stand qualitatively the effect of this term, we neglect 
momentarily other terms and approximate H p h by a con- 
stant ( we recall that a static bubble n p /i(q, 0) is a con- 
stant for q < 2fcp in 2D). After simple manipulations, we 
then obtain 



1 



1 



i + u 2 u 2 ph 2 1 - un ph 



u 1 

"2 1- UU ph (K - P 



■&a8 ■ CTp-y- 



(A5) 



We see that the term 1/(1 + UU p h) describing the inter- 
action via soft bosons in the charge channel is no longer 
there. Although disappearance of this term well may 
be an artifact of the approximation [Eq. (|A5[) is, strictly 
speaking valid only to order U 3 ], this is still clearly a 
warning sign for the whole approach, as it shows that 
additional terms with particle-hole bubbles, not included 
into RPA-type analysis may be relevant. Note also that 
there is a "leakage" of the 1/(1 — UH p h) term from the 
spin channel into the charge channel, i.e interaction me- 
diated by a soft boson in the spin channel induces the 
same interaction in the charge channel. This effect was 
considered in detail in Rcf. \2(1 



APPENDIX B: ALTERNATIVE AVERAGING 
PROCEDURES FOR T n 

For completeness, we present here the results for Zr, 
Z, and to* jm obtained within alternative approaches for 
summing up the diagrammatic series (|3.30p for T n . We 
remind that the uncertainty is related to the presence of 



the factors d] fL d p (d 



■2(«-l) , j2(n-l) 



) in the diagrammatic 



series for T^ 7 ,$(k, p). In the main part of text, we used 
the approximation in which the angular dependence of 
^a/3;-yS was averaged over the FS. An alternative is to 
keep only the d^d p term at each order and neglect other 
(non-d— wave) terms. A simple exercise in trigonometry 
shows that this amounts to replacing 



j , / ,2(n-l) 

«kctp(a k 



(2n - l)rf k d p (Bl) 



at each order of the perturbation theory. Summing up 



the series for Y apnS 
valid, i.e., 4 = T 



n{n} (k,p), we find that Eq. $£5Q is still 
ap-yS /Zr, where Zr is a constant, 



but now the expression for Zr is different: 



1 

— = 1 

Z T 



OO 

E 

n=l 



Z 2 m* 



2to 



/ -i _ Z 2 m'x \ 
y 1 2m J 



2 1 



2m t 



(B2) 



Using this form and isotropic FL formulas for to*/to and 
we obtain after simple algebra that in the critical FL 
regime, where A is large, m* jm = 1/Z = A/2 and Zr ~ 

l/2A«Z/2. 

We see that, as in Eq. (|4.18[) . ^r| A>1 ~ Z, the only 
difference between (|4.18[) and the present case being the 
numerical prefactor. This proves our point that both 
procedures for summing up the ladder series for T n yield 
physically equivalent results. 

We also obtain very similar results by ne- 
glecting the d— wave structure of the intermedi- 
ate vertices in the ladder scries, i.e., replacing 



,2(n-l) 



)/2 by (1/2)' 



Z 2 m*\ 
4m 



In this situa- 
). Combining this 



tion, we obtain 1/Zr = 1/(1 

with the equations for Z and m* /to and solving the full 
set, we obtain to* jm = 1/Z = A/4, and ^r| A>1 = Z for 
A> 1. 



APPENDIX C: LOOP EXPANSION FOR E 

a. Relation between vertex renormahzation and 
single-particle residue 

We begin by discussing a subtle point in the rela- 
tion between the effective fermion-boson interaction D = 
9fbX an d ZF n . We found that in the critical FL regime, 
the vertex renormahzation constant Zr is related to the 
single-particle residue Z as Zr = CZ, where C ~ 1 
depends on which of the approximate averaging proce- 
dures, discussed in Sec. IIII CI and Appendix [B] is em- 
ployed. One of such procedures yields C = 1, whereas 
the other two yield C ^ 1. The uncertainty is related 
to the fact that we used anisotropic d— wave form of 
the interaction between fermions, yet approximated the 
quasiparticle residue Z and the effective mass m* by 
constants independent of the position on the FS. Re- 
gardless of a particular value of C, proportionality be- 



f2,RPA 



is 



tween Z r and Z implies that ZT U = (Z/Z r )T : 
proportional to r^' 11 ^. It is then instructive to ver- 
ify whether the exact Pitaevski-Landau formula for £ is 
reproduced if the diagrammatic expansion for the ver- 
tex T n in terms of r n ' RPA is transformed into an ex- 
pansion for the self-energy by contracting a pair of ex- 
ternal legs of a vertex. The first term in the series is 
r n,RPA , and hence the "bare" self-energy is the one with 
r o,RPA i nstea d f D. This would be the final result if 
Z and Zr were equal in the critical FL (in that case, 

x = zr n = {z/z r )T n > RPA - r n > RPA ). 
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Whether Zy is equal to Z or just proportional to Z 
does not affect our main result that all Landau parame- 
ters diverge upon approaching a nematic QCP. It turns 
out, however, that only an exact equality Zp = Z is con- 
sistent with the loop-wise expansion for the self-energy. 
This is not surprising because an equality Zy = Z is the 
result of a particular averaging procedure in which the 
<i-wavc vertices arising from intermediate states were re- 
placed by constants when evaluating the ladder series for 
r n . This is consistent with replacing Z and to* by angle- 
independent constants. In the other two procedures, we 
either averaged the angle-dependent interaction over the 
FS at each order or extracted the d— wave component 
from the interaction. In both cases, there are extra com- 
binatorial factors associated with the d— wave nature of 
the interaction. 

To see that only an equality D = r n,RPA is consistent 
with diagrammatics, we consider explicitly diagrams for 
S. To second-loop order, there are two diagrams, b and 
c in Fig. [7] Diagram b is a part of the self-consistent 
renormalization of the Green's function in the one-loop 
diagram. Diagram c is a vertex correction to the one-loop 
self-energy In two subsequent Sections, we will compute 
these two diagrams and show that they are indeed small, 
at least as along as akp ^> 1. 



b. Green's function renormalization in the one-loop 
self-energy 



We start with diagram b in Fig. [7] 



S 2L ~ / d 3 Pd 3 Q [G 2 (P)] n G(Q)D(K - P)D(Q - P). 

(CI) 



Integrating over Q first, we obtain 

4l ~ J d 3 P[G 2 (P)] u D(K-P)XW(P), (C2) 

where Sii(P) is one loop-self-energy with D = r n ' RPA . 
The leading term in £ix(P) is Xuj. This term, however, 
vanishes upon integration over P. The most straight- 
forward way to see this is to perform momentum inte- 
gration first. Replacing J d 2 p in Eq. (|C2[) by J de p and 
integrating over e p from — A to A, where A ~ E-p, we find 
after simple algebra that, because of the double pole in 
G 2 (P), the integral comes either from vanishingly small 
e p or from \e p \ ~ A. The high-energy contribution to 
E^ 2 } is of order 



A / dbJr, 



det 



c *\2 



n 2A lu 2 



The contribution from vanishingly small e p 



Z 2 m* 



-A 



m 



ddJpUJp 



del 



id sgne p ) 2 



(C3) 



(C4) 



comes from the branch cut at e p = 0, caused by sgnej 



term in Eq. (|C4[) . Evaluating the integral, we find that 
it vanishes in the limit of 8 —> 0: 



Z 2 to* 



-A 



dWqUJq - 



(C5) 



A non-zero contribution to X^z? comes only from the ek 
dependent term in Sil, given by Eq. 1)4. 20j) : 



J 2L 



Z 2 Ar 



m 



dujpd p 



(w - e* +iS e } 



;D(K~P) 



ZT/^m* 



duj p d 2 pG(P)D{K - P) ~ TaEil. 



(C6) 



r 



As we have shown in Sec. IIV CI Ta ~ l/akp <C 1 i.e., 

is indeed small compared to Ei£. The same con- 
clusion holds for higher-order diagrams of the same type. 



c. Vertex corrections to the one-loop self-energy. 

Next, we consider the vertex correction to the one-loop 
self-energy given by diagram c in Fig. [JJ Integrating over 
P, we obtain 



4l(K) = -J -^r A (Q;K)D(Q)G(K + 



Q). (C7) 



In Sec. MI C 3| we computed the three-leg vertex 
T/\(Q;K) in the limit of il/q — > 0. In this limit, 
Ta ~ 1/akp <C 1 coincides with the effective Migdal 
parameter of our theory. To calculate the self-energy, 
however, we need to know Fa(Q;^) at arbitrary £l/q. 
When evaluating the dynamic part of the vertex, it is 
important to account for finite curvature of the Fermi 
surface, i.e., to keep a quadratic term in the expansion 
of the dispersion 



t k P +q 



v fQ\\ 



2m* 



(C8) 



where gii and q± are the components of q along and trans- 
verse to kp, correspondingly. In general, Fa(Q; K) is an 
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a) 



b) 




K K+Q 2 K+ Ql+Q 2 K+Qj 



K 



FIG. 7: Diagrams for the self-energy at the one-loop [a)] and 
two- loop orders [&) and c)]. The wavy line is the RPA inter- 
action, Eq. ([3T23]) . 



involved function of its arguments; however, a useful es- 
timate can be written as 



r A (q,tt;k F ,w fe = 0) = r A +A 



Z 2 m* 



101 



m max{|Q|,|e£ F+q |}' 
(C9) 

where the first term is the vertex in the limit fl/q — > 0. 
As it is always the case in Migdal-Eliashberg-type the- 
ories, the vertex depends crucially on the ratio of en- 
ergies entering the dynamic term. Substituting the dy- 
namic term into the self-energy, one can readily see that 
f2 is of order of the external frequency uJk, q» is of order 
uj /v F , and q± is of order of the inverse correlation length 



£ 1 = (1 + gca) 1 ^ 2 ' l a i so that a typical value of |e£ p+q | 
is the largest of the two energies- \uj^\ and wfli where 



1 



3/2 



m*£ 2 

Therefore, the vertex can be estimated as 



r A (q,fi;k F ,^ fe = o)~r 



mi 



A 



max{|fi|, wfl} 



(CIO) 



(Cll) 



The physical meaning of ujfl is that it is an upper bound- 
ary for a FL behavior near a QCP. For energies higher 
than ujpL? Landau damping becomes the dominant term 
in r n,RPA and, consequently, Z and m* become energy 
dependent. The FL theory is therefore valid only at ener- 
gies smaller than wfl- We see from Eq. (|C11|) that vertex 
corrections are irrelevant for La <C 1 and \lo\ <C ujfl- 

Notice that the smallness of dynamic term in L A is 
due to the presence of the quadratic term q^_/(2m*) in 
the fermionic dispersion, Eq. (|C8[) . This term reflects fi- 
nite curvature of the FS in larger than one dimensions. 
Without such a term, vertex correction would be of order 
one for any frequency. Importance of the the FS curva- 
ture for vertex corrections in 2D has been discussed in 

Refs. nMnnomm. 

To verify this reasoning, we computed diagram c in 
Fig. [7] explicitly. There are two contributions to the self- 
energy: the first one is linear in uj and gives a correction 
to mass renormalization; the second one is quadratic in 
uj and gives a correction to damping. Although the uj 2 
contribution is formally smaller than the linear one, the 
coefficient of the uj 2 term diverges in the non-FL regime. 
The linear contribution can be extracted by keeping only 
L A in ([U9|l and substituting it into Eq. (|CT)) . This gives 



£<f } (k F ,u,) 



La^u 



LaAw. 



(C12) 



The uj 2 contribution is obtained from the entire expres- 
sion for the self-energy 



£^' 2} (k F ,w) - Z 3 J d 2 q 1 d 2 q 2 dn 1 dtt 2 



1 n. rpa 



( qil n 1 )r^ A (q 2 ,n 2 ) 



i {uj + f2i) 



k F +qi+q2 



i(oj + n 2 ) 



-k F +q 2 

(C13) 



where e^F+q ^ s gi ven by Eq. (|C8[) . Integrating over <j»|| 
first, we see that the region of integration over the in- 
ternal frequencies, f^i and Q 2 , is bounded by external 
uj. This immediately shows that the double integral 
J J dVlidQ, 2 contributes a factor of uj 2 . Rescaling vari- 



ables as x = and y = q 2 ±£, we obtain 



v {2,c}/, N -xl^l 3 [°° dx f°° 



dy 



1 



+ y 2 (3 2 + x 2 y 2 ' 
(C14) 

where j3 = \uj\/uj-p\ J . For (3 <C 1, the double integral 
in Eq. (|C14|) behaves as ln|/3|/|/3|. Collecting the two 
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contributions, the final result for diagram c reads 
V { C } / vp/a , . w \ 

We see that T,^ is indeed parametrically smaller than 
Sil in the FL regime, i.e., vertex corrections are also 
irrelevant for the fermionic self-energy. 

To conclude this Section, we compare our results with 
the canonical Migdal-Eliashberg theory of the electron- 
phonon interaction. In that theory, the vertex correction 
to the self-energy has a similar form, except for the de- 
nominators of both terms in Eq. (|C15[) contain the same 
energy scale: the Fermi energy. In this case, the sec- 
ond term is always smaller than the first one as long as 
lu is smaller than the typical phonon frequency. In our 
case, however, the effective "Fermi energies" in the static 
and dynamic terms are different and, what is most im- 



portant, they behave differently as a QCP is approached: 
whereas E-p remains finite, the upper boundary of the FL 
behavior, wfl, vanishes as 1/A 3 . In the non-FL regime, 
the one-loop self-energy behaves as Six ~ u)^ 3 u> 2 / 3 for 
uj -C Woi where ojq = E-p/{akp) i . This behavior can be 
interpreted as resulting from the energy dependences of 
the mass renormalization coefficient X(uj) ~ (luq/uj) 1 ^ 3 
and of the quasiparticle residue Z(u>) ~ 1/A(cj). Re- 
writing wfl as loq/X 3 with w-dependent A, we see that 
~ id. Therefore, the second term in Eq. (|C15[) is of 
order unity, i.e., vertex corrections are, in general, im- 
portant in the non-FL region of a QCP and also for the 
uj 2 In a; term in the FL regime. 3 - Some of the higher- loop 
diagrams can be regularized by expanding the theory to 
the large N cas o 25 ' 28 i 33 ; however, there are still diagrams 
that arc not small in the large N approximation. 38 ! 39 
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